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End  fltresa  calculation*  on  Elaatic  cylinders 
by  P.J.D.  Mayas  and  D.A.  spanca 

ABSTRACT 

For  a  semi- infinite  circular  alaatic  cylindar  s>0,  r<l  daformad 
solaly  by  a  distribution  of  stress  and  displacamants  on  its  flat 
and  z»0,  the  Love  strasa  function  can  ba  axpandad  in  a  sarias  of 
eiganf unctions  of  known  form.  For  problems  in  which  mixed  stress 
and  displacements  boundary  conditions  ara  prescribed  on  s*o  the 
coefficients  appearing  in  the  expansion  can  be  determined  in  an 
explicit  form  via  seta  of  biorthogonal  functions.  When  normal  and 
shear  stresses  are  prescribed  on  z*0  no  such  closed  expressions 
for  the  coefficients  exist  and  approximate  methods  usually  lead  to 
infinite  systems  of  linear  equations  which  are  solved  by 
truncation,  stability  of  solution  as  the  order  of  truncation  is 
increased  can  only  be  guaranteed  theoretically  when  the  infinite 
matrix  is  diagonally  dominated,  and  this  is  not  the  case  for 
existing  methods.  A  Galerkin  method  has  been  developed  using 
weighting  functions  chosen  so  as  to  optimise  the  diagonal 
dominance  of  the  infinite  matrix,  and  numerical  results  show  that 
although  the  resulting  matrix  is  not  completely  diagonally 
dominated,  the  resulting  coefficients  show  an  improvement  in 
stability,  and  accurate  solutions  can  be  obtained  using  smaller 
matrices  thus  producing  a  much  more  efficient  method  of  solution. 
Calculations  are  presented  numerically  and  graphically  for 
representitive  distributions  for  three  classes  of  datai- 

(i)  Smooth  continuous  data 

(ii)  Smooth  data  violating  compatibility  at  r>l 

(iii)  Data  containing  discontinuities 
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End  stress  calculation*  on  Elastic  Cylinders 

i 

by  P . J . D .  Mayas  and  D . A .  Spanca 


The  Love  stress  function  $(r,s)  in  an  alastic  cylindar  *>o, 
r<l  subjected  to  homogeneous  boundary  conditions  on  the  curved 
boundary  r» l  can  be  expressed  as  an  eigenfunction  expansion  of  the 
form 


»'  X"%(r»Xn) 


(1.1) 


where  X  is  an  eigenvalue  determined  from  the  conditions  on  r*l . 
n 

For  the  case  of  a  traction- free  curved  face,  Xn  is  a  root  of 


X2{jJ(X)  +  jJ(X)}  *  2(l-V)jJ(X) 


(1.2) 


Little  and  Childs  [1967]  have  given  a  construction  for  determining 

the  coefficients  c  in  the  expansion  (1.1)  for  cases  in  which  the 

n 

data  on  the  flat  end  z=0  takes  the  form  of  prescribed  values  of 
either  of  the  pairs 


or 


ur 

°r.  ■nd 


(1.3) 


For  these  "canonical**  problems  the  {c^}  are  found  explicitly  as 
quadratures  of  the  data  with  appropriate  biorthogonal  functions 
derived  from  the  +(r>Xn). 
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in  the  present  report  we  consider  the  problem  of  determining 

the  coefficients  when  <J  and  o„_  are  prescribed,  it  is  known  that 

zz  rz 

no  explicit  solution  exists  for  this  case,  and  the  (cn)  must  be 
found  by  approximate  methods  leading  in  general  to  infinite 
matrices  which  can  only  be  inverted  in  truncated  form. 

This  problem  has  already  been  studied  extensively  for  the 
elastic  strip,  x>0,  | y | <1 .  Spence  [1978]  introduced  a  set  of 

weighting  functions  derived  from  members  of  the  family  of 

'  *  V. 

biorthogonal  functions,  which  in  the  case  of  the  traction 
problem  for  the  strip,  namely 

a  .a  defined  on  x= o 
xx*  xy 

lead  to  a  diagonally  dominated  system  of  equations 


TK  c  *  d 
^  mn  n  m 
n 


(1.4) 


where  K  *  I  -  G,  with  the  row  sum  norm  IGI  <  1.  For  such  a  system, 
IN) 

the  solution  c  ’  say,  of  the  truncated  system 


n 


<N)  (M) 
mn  n 


d<N> 

m 


(15) 


is  known  to  converge  to  the  solution  oi  the  full  system  as  N  -  «», 
and  this  was  borne  out  for  the  cases  tested,  in  which  it  was  found 
that  changing  the  order  of  truncation  N  did  not  lead  to  changes  in 
the  coefficients.  This  was  not  found  to  be  the  case  with  other 
published  methods  that  were  tested. 


a  - 


2.  The  Nwr  formulation 


The  construction  given  by  Little  and  Childs  [1967]  for 
obtaining  biorthogonal  functions  for  the  two  canonical  end 
problems  for  the  elastic  cylinder,  thus  enabling  them  to  obtain 
the  coefficients  appearing  in  (1.1)  explicitly,  has  not  proved  to 
be  the  most  suitable  for  the  present  studies.  The  main  disadvatage 
is  that  for  the  stress  problem  it  is  not  possible  to  "optimise” 
the  weighting  functions ,  thus  improving  the  diagonal  dominance  of 
the  infinite  matrix  arising  in  this  problem.  Consequently  we 
choose  a  different  but  equivalent  set  of  four  stress-  and 
displacement- related  variables  which  will  be  prescribed  on  x*0. 

In  terms  of  the  biharmonic  "Love"  stress  function  (Love 
[1927],  Art. 188)  the  stresses  and  displacements  are  given  by 


„  -  fee  «-»"*• -s»  «•*•»» 


62*1 


zz  =  §j{(a-v)v2*  -  HI)  °eo  *  '  ?!?)  (2.3,4) 


32A  I  32$ 

2Mur  =  *  2(1-V)V  »  - 


(2.5,6) 


where  v  is  Poisson's  ratio  and 


_2  13  3  a*  _  .  a* 

v  x  +  IJI  *  B  +  TZT 


(2.7) 


if  the  cylinder  is  subjected  to  stress-free  side  conditions  on  r«i 
and  a  self- equilibrating  distribution  of  stresses  and 
displacements  on  «»o  then  ®  may  be  expanded  as  an  eigenfunction 
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expansion 


4>(r ,  z)  =  Ecm# (r  t  XB) e"  Xm*  (2.8) 

m 

where  X_  is  a  root  of 

HI 

Xa {j*  (X)  +  jJ(X)}  -  2(l-»)jJ(X)  »  0.  (2.9) 

♦(r»X)  *  |2<l-v)Ji(X)+Xai(X)]j#(Xr)  +  XJ4 (X)rJj(Xr)  (2.10) 


and  the  correct  interpretation  of  the  summation  (2.8)  is  obtained 
by  numbering  the  roots  of  (2.9)  in  the  right  half-plane  so  that 
1  s  x  [see  figure  1]  and  writing  the  expansion  more  precisely 


4>(r,z) 


«=+«* 

E'cm4<r» 

«p-« 


x  >< 


■  x_* 


(2.11) 


where  the  prime  means  that  the  term  with  m*0  does  not  appear  in 
the  summation .  This  implies  that  the  normal  stress  distribution  is 
equilibrated,  i.e.  (lra _ (r.O)dr  *  0. 

The  present  choice  of  prescribed  functions  together  with 
their  expansions  in  terms  of  the  "derived”  functions  +^a> (r)  are 
given  by 


(r) 

a°«s/dr 

♦m”  <r> 

*'a> (r) 

cr. 

(r> 

<*> 

m 

<r> 

*'4' (r) 

>■0 

C*  <r> 

This  can  be  seen  to  be  equivalent  to  prescribing  the  unmodified 

stresses  and  displacements  as  in  Little  and  Childs  -  for  example, 

if  a  and  u  are  known  on  z=0 ,  then  so  are  f ' n  and  f19*  as 
zz  r 

defined  above. 

In  terms  of  $(r;X)  the  derived  functions  ♦'an  are  given  by 

ut 


«*  ■  <^h*2*  4  -<■£ 


.  (  3  >  .  . 
♦m  <r> 


-Xm{-(X-2W)|?B%  +  2VX*.ff} 


C’<r>  *  xi.fi} 


m  drJ 


(2.13) 

(2.14) 

(2.15) 

(2.16) 


and  explicit  expressions  for  these  functions  in  terms  of  Bessel 
functions  are 


C<»  ■  Ck»wrw>  4 


♦«”<*>  *  kJ{-WrVx«r>  4  WW>) 


(2.17) 
(2. IB) 


<*'  <*> 


.  «*)  .  . 
♦m  <r> 


X^X*VXm>rVXmr>  +  <Xm>+X»J,  <X»>  JJ1  <X»r>  ) 

(2.19) 

-a(l-n>)»’j1(X„>Jl(X1,r) 


(2.20) 


4>(rjX)  is  a  solution  of  the  reduced  biharmonic  equation 


ri  5L_r!_ 

lr  dr  dr 


(3.1) 


and  as  in  Spence  [1978]  this  equation  may  be  expressed  as  a  matrix 
differential  equation  in  either  ♦,l>  and  or  ♦*a>  and 

RI  Rl  In  In 

For  the  (1, 3) -canonical  problem  the  matrix  equation 


-(2+V)[Ba-ij 


«1» 

m 

.  (3) 


a-...)  ^ 


(3.2) 


can  readily  be  shown  using  (3.13,15)  to  reduce  to 


(3.3) 


The  condition  o  =  o  on  r*l  may  be  written  in  terms  of  + 
rz  i 

♦  ‘3’  as 


en 


and 


l*'11  (1)  ~  ♦  ***  (1)  -  0  (3.4) 

m  Rl 

The  corresponding  boundary  condition  for  ofr  is 
(l-v)D#Jk1’ (1)  +  2D^3,(i)  +  <i+v)^3’(l)  *  0 


(3.5) 


where  D  •  d/dr.  The  derivatives  of  4' 11  and  4  131  contain  the 

in  hi 

fourth  derivatives  of  4,  and  in  obtaining  (3.5)  it  has  been 
necessary  to  use  the  reduced  biharmonic  equation  (3.1)  evaluated 
at  r=l  to  express  the  orr  condition  in  the  required  form. 

As  in  Spence  [1978]  the  function  1 1  and  f^3  *  which  are 
biorthogonal  to  4 1 1 *  and  4 ' 3 1  are  obtained  as  the  eigenfunctions 

HI  HI 

of  the  differential  operator  adjoint  to  (3.3)  which  are 
constructed as  follows  1 - 

Using  the  differential  equation  (3.2)  we  may  write 


■  <f  <»’>*:*' K" 

*  ((-(•S.lC’  ♦ 


"■o  ( <*>  ■  /:*•«“  ] 

W“ 

(a+.)  -  <>> 


we  may  now  integrate  twice  by  parts,  transferring  the 

B3  derivatives  onto  the  f ,an  and  introducing  boundary  conditions 

n 

at  r=li- 


1 

1 

i 

i 

j 


The  construction  of  biorthogonal  functions  for  the 
(2,4)-problem  is  a  modification  of  the  work  of  Klemm  [1970],  who 
treated  the  full  non- axi symmetric  end  loading  problem.  Putting 
0eo,  a/ae-0  in  his  construction  gives  the  biorthogonality  given 
here.  However,  his  construction  for  the  (l, 3) -problem  does  not 
lead  to  a  pure  biorthogonality  from  which  the  coefficients  can  be 
determined  explicitly,  and  the  construction  described  below  is 
new. 
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-  (O’-I.K")-  <**'>(<" l''-hK3,> 


(l+v)<kM-Xn)<4m  tn  +  ♦.,*»>  “ 

+  |(i+v)4^1’ (i){-v^11  <i)  +  a(i+»)^*'a)  +  v(a>»)^*’  (i)} 

♦  (Df-*"' <i>  -  •»£**<*>)  (3-7) 

Thus  if  3> ■•tisfiss  ths  adjoint  boundary  conditions 

(1)  +  Vf^’U)  »  0  (3.8) 

-v,“*  (X)  +  aa+v)^*’ <i)  +  vta+vjf^’d)  ■  o,  0.9) 


or  mors  compactly 


9 


(3.10) 


K"  (1>  +  U> 


D^3’  (1)  -  (1) 


=  0, 


(3.11) 


we  find 


wuc  +  co  •  ° 


(3.12) 


and  hence 


<cc  -  co  - °  — • 


(3.13) 


Exactly  the  same  conatruction  may  be  used  for  the 
(2, 4) -canonical  problem.  This  time  the  required  matrix 
differential  equation  is 


-(1+v)  [b* 

-J,l  ["'-Ill 

|C  1 

.  2 

Cl 

o 

-(l+v) (oa-72J 

lei 

-  <!«■>*» 

cl 

(3.14) 


with  corresponding  boundary  conditions 


(3. IS) 


(l+V)0^21  (1) 


DC  <*>  +  <*'  (1) 


(3. 16) 


and  the  adjoint  equation  and  boundary  conditions  are 


(« 


U+v)[B3-i,]  o  |K3‘ 
S,1  -(itv)fB3-ia]  U-v 


d+v)x; 


rn** 

*«4» 


(3.17) 
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(3.18) 


(l+v)©*"'  (l)  =  D^21  (l>  +  v^a’  <1) 

(1)  =  O  (3.19) 


resulting  in  the  biorthogonality 


<♦, 


f  21  f21  +  .  141  (41 

m  *n  *m  ”n 


>* 


for  m»*n. 


(3.20) 


in  terms  of  the  Bessel  functions  the  two  biorthogonal  vectors  are 
given  by 


tn  (r) 

=  A 

^.C  31  .  . 

[*n  (r) 

n 

0-21) 


<2'  <*> 

*  B 

2(i+v)J1<Xn)<»4anr) 

*n'  <r> 

n 

WV'VV^nVVVV^ 

where 

A  =  _ l _ 

n  2 2  <*„>*(*„> 

B  =  ...  -i-  — 

n  2(i+v)lnJia<Xn)P(ln) 


(3.22) 


(3.23) 

(3.24) 


*<x„>  *  ~knJllk  n>  +  a<1~v>  XnW  VXn>  “  ad-ajjja^  <»•»> 


and  the  normalising  factor  P ( ln)  has  been  introduced  so  that 


/xiu.m  .  .  (* i c» i \  . 

(♦m  *n  +  ♦«  *n  >  “ 


(3.26) 


-  10  - 


+  ♦ 


(41 

m 


mn 


(3.27) 


it  is  interesting  to  note  that  as  in  Spence  [1978]  this 
formulation  exhibits  what  might  be  called  a  "self-biorthogonality" 
where 


fr'ul 

wn 

A 

n 

(S. 

lFn 

(l+v)Xn« 

(l-V)^n  +  2^s* 

and 


rn 

-  B 

n 

c* 

(4’ 

lfn 

2  (1+V) X  » 

HI 

2  (  !+ V)  4j|2  * 

(3.28) 


(3.29) 


By  contrast,  in  the  formulation  of  Little  and  Childs,  the 

(1. 3)  -biorthogonal  functions  are  given  in  terms  of  the 

(2. 4)  -derived  functions  and  vice  versa. 


11 


4.  Optimal  Weighting  Function* 


In  this  section  we  consider  the  stress  problem  in  which 


and 


(4.1) 


^r,).=  0  "  ',a‘<r> 


are  prescribed  functions  of  r.  This  does  not  fall  into  the  class 
of  canonical  end  problems  categorised  in  Section  1.  As  was  done 
for  the  strip  problem,  we  now  seek  weighting  functions  of  the  form 


X^1’  *  +  "♦m** 

*1**  -  c*iC  +  °*W4‘ 


(4.  a) 
(4.3) 


where  A,  B,  c  and  D  are  constants  to  be  determined.  [The  choice 
A  *  -2v,  B  *  (l-v) ,  C  ■  0,  D  ■  -(14V)  would  produce  the 
biorthogonal  functions  ,  f™1  defined  in  section  3,  but  as 
will  be  seen  these  are  not  optistal  for  the  non- canonical  problem] . 

An  infinite  set  of  linear  equations  for  the  coefficients 
cr  in  the  derived  expansions 


(!) 


121 


Ecn*na> 


(4.4) 

(4.5) 


is  obtained  by  combining  the  scalar  products  of  (4.4)  with  xJ^U 
and  (4.5)  with  x^a>  for  each  n.  This  yields  the  set 


-  la  - 


VK  c  =  d 
mn  n  tn 


where 


(46) 


(4.7) 


end 


t 

we  now  chooee  the  constants  A,  B,  c  and  D  so  as  to  make  the 
off-diagonal  elements  of  the  matrix  A  as  small  as  possible  in 
absolute  value  compared  with  the  diagonal  elements.  Por  this 
purpose  the  scalar  products 


<.<“<")■  (CO-  (CO-WY)  <*••> 


have  been  calculated  and  are  listed  in  Appendix  A.  The  expressions 

are  cumbersome,  but  the  salient  feature  is  that  the  first  three 

contain  the  factor  (Xa-X*)“*.  As  was  noted  by  Spence  for  the  strip 

■1  n 

problem,  the  presence  of  any  negative  power  of  (XM-Xn)  in  the 

matrix  A  leads  to  divergent  row  sum  norms.  The  four  constants  A, 
mn 

B,  c  and  D  provide  just  sufficient  freedom  to  suppress  all  such 
factors  in  the  denominator . 

The  procedure  for  determining  the  optimal  choice  for  the 
constants  A,B,c  and  o  given  the  choice  of  weighting  functions 
(4.2,3)  involves  taking  the  matrix  elements  (4.7)  with  and 

X,a*  given  by  (4.2,3),  and  dividing  out  the  unwanted  factors 
giving  three  equations  for  the  four  constants. 


-  13  - 


Using  th«  quadratures  given  in  Appendix  B  the  general  matrix 


element  A  is 


*«,  *  *«“C>  -  '«"<")  *  ♦  WO 


1,  A- <**»»>  X.X„<X»VX«>'’1 <X»>-X»VX«>  Vx»>  >1 

1  n  n 


♦  O  d[  <*+■»'> > 

m  n 


+  “<x+»>xixiVx«>VXr.>  -  oxixn<x.J,<x«>,.<x„>+xnJ.<x«>',,<x„>> 


♦  D  <x«>  xixn«V.  <x-> J.  <  V-VW  J.  <V>1 
+  (VV**1  »  <M>  x^  W  V  V  V  *  W  V  V  V  > 

~x:*x.vx.>vv-x»vx.>vv>)} 


♦  <<x-')xixi,;(x.>,!<x»>[i-irrs(-»<«|’>x»x„-o<x«)x;) 

■  n 

♦  TT^. l*x,xn < X> x^>  -»x,x, < x> x^> ■« < xi« >n>  xil 

we  shall  try  to  eliminate  the  (X^- X*)‘  1  from  the  dominant  term, 
which  over  a  common  denominator  can  be  written 

«xixivx.)j.  <xn»  ( f-  «!-»<x+*>  >  xi*»  <»+«>  xix, 

ww 

-  (C4D  (Its)  )  *-  (»+»*)  1*X  *♦  (A+BV)  A*1  *J  am)*$an> 
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+  [-  (A+BV)X^4(C-D(l+V))X*Xn 

+  ( 3A-2B+BV) X  *  X  *4 (C4D (14V) ) X  * X  * J  J# ( Xm)  ( Xft) J 

The  condition  that  both  factors  multiplying  the  Bessel  functions 

have  a  factor  X  -X„  is  the  same,  namely 
in  n 

A  -  B  4  C  *  0,  (4.10) 


and  if  this  condition  is  satisfied  the  dominant  term  becomes 

<Xm)  (Xn}  (Ji  (Xm)  J*  <Xn>  f  (C-D<l+V>  >  (2A42BV4C4D  (14V)  )  X^ 


<wW 


- (2A42BV) X *X *4 (A4BV) X^X *4 (A4BV) X*] 

+  XmJs  <Xm>Ji (XnJ  f  <A+Bw>  x«+ (3A-3B4BV4C4D (14V) ) XmXn4 (C4D(14V) ) X*] J 


Again  the  condition  that  both  terms  inside  the  square  brackets 

have  the  factor  X  -X„  is  the  samei 

m  n 


A  4  BV  4  O ( 14V)  «  0 


(4.11) 


giving  a  dominant  term 

-♦‘M*.  's.>w  (•»,  f  (»♦»«««  »»»» 

<w*<w 

- (2A48BV) X^X*- (A4BV) X*] 


IB 


+  XmJo  (Xm>  Ji  <V  f  <A+BV>  X»“  <C+D(l+v>>  xj  J 


The  last  relation  suppressing  all  factors  1  in  the 

denominator  is 


A  +  BV  -  C  -  D(l+V)  =  0 


(4.12) 


giving  as  the  dominant  term 


«<~->*::vvvv{<3*’w^vvvv^vvvv) 


<X»*X»>’ 


The  three  equations  (4.10,11,12)  lead  to  the  values 


A  *  -(1-2V),  B  a  -3»  C  a  -2(l+v)»  D  «  1. 


(4.13) 


The  resulting  weighting  functions  are  thus 

(4.14) 

C  -  2<l+»x«lx,.Wr'W>  -  n.i'.i'W'.ii'.'Vl 

(4.15) 

and  the  matrix  elements  are 


4  <*♦'>  <XJ  VV  (<,xi*sV»*x^  vv  vv^ivv  Vx»> ) 


(4.16) 


\m  *  4  W*.’)  <<17> 

In  order  to  see  why  this  choice  of  coefficients  should  give 
rise  to  a  more  stable  matrix,  it  is  of  interest  to  determine  the 
asymptotic  form  of  the  matrix  elements.  These  are  determined  by 
using  the  asymptotic  form  for  the  eigenvalues 

X  “  pi«  +  4ilog(4ms)  (4.18) 

m  * 


so  that  we  simply  replace  X^  by  ms  to  first  order  in  the  matrix 

elements  (4.16,17).  Using  the  asymptotic  forms  of  the  Bessel 

function*  given,  for  example,  in  Abramovitch  and  Stegun,  it  can  be 

shown  that  the  Bessel  functions  of  the  eigenvalues,  namely  J  (X  ) 

i  ® 

and  J  (X  )  ,  have  the  asyinptotic  form 
i  m 


w 


(-1) 


illil 

Vs 


Vx»>  '  < 


-l)"*1  lip LI 


(4.9) 


Thus  the  matrix  elements  have  the  asymptotic  form 


Atnn  * 

for  X  in  the  first  quadrant,  and 
in 

Amn  ‘  . 

for  X  in  the  fourth  quadrant,  with 

m 
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r’ 

A _  “  8{1+V) «*®* 

vran 

if  the  factors  (X  -X  )  1  had  not  been  eliminated,  then  the  row 

ri  n 

sums  would  grow  with  m,  for  exactly  the  same  reasons  as  given  in 
Spence  [1978]  for  the  strip  problem.  , 
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s. 


In  order  to  test  the  optimal  weighting  functions  derived  in 
section  4  and  compare  them  with  unmodified  biorthogonal  weighting 
functions  the  following  sample  stress  distributions  were 
considered 


19  - 


In  order  that  the  stresses  should  decay  as  z-«*>  the  normal 
stress  must  be  self-equilibra-ced.  All  the  distributions  tested 
satisfy  this  condition.  For  non- self-equilibrated  distributions  a 
simple  polynomial  term  can  be  added  to  a  stress  function 
representing  am  equilibrated  distribution.  in  addition  to  this 
condition  on  the  normal  stress,  the  shear  stress,  as  well  as 
varnishing  at  the  origin,  must  also  be  zero  at  r=l  if  the  end 
distribution  is  to  be  compatible  with  the  zero  stress  condition  on 
r=  1 . 

The  first  two  cases  satisfy  the  conditions  of  equilibration 
and  compatibility  and  are  continuous.  Cases  3  and  4  have  simple 
jump  discontinuities  in  the  first  derivative  of  the  prescribed 
stresses.  Case  S  is  equilibrated  but  discontinuous,  and  case  6  is 
incompatible  with  the  side  conditions  on  the  shear  stress. 

It  is  only  possible  to  find  closed  forms  for  integrals  of  the 

form 

when  k  is  even  for  the  integrals  and  k  is  odd  for  the  J ^ 
integrals.  Therefore  the  prescribed  normal  stress  distribution  may 
only  contain  even  powers  of  r  and  the  shear  stress  distribution 
odd  powers,  if  the  right-hand  sides  of  the  truncated  systems  are 
to  be  evaluated  in  closed  form.  Using  integration  by  parts  the 
other  integrals  may  be  reduced  to 

fw>dt 

*  0 

-  ao  - 


which  could  be  evaluated  numerically.  However  the  real  and 
imaginary  parts  of  J#(Xnt)  become  more  oscillatory  as  n  increases, 
which  presents  problems  for  library  integration  subroutines . 
Although  a  general  program  for  solving  the  end  stress  problem 
would  need  to  include  the  possibility  of  general  polynomial  stress 
distributions,  for  the  purposes  of  this  report  it  was  decided  that 
sufficient  test  could  be  devised  with  the  above  restrictions. 

Three  salient  features  of  the  numerical  results  presented  in 
appendices  C  and  D  are  worthy  of  note,  showing  the  advantages 
of fered  by  Optimal  Weighting  functions .  These  are 

(i)  The  improvement  in  diagonal  dominance  of  the  truncated 
matrices . 

(ii)  The  increase  in  stability  of  the  earlier  coefficients  as 
the  order  of  truncation  is  increased. 

(iii)  improved  convergence  to  the  data  for  various  orders  of 
truncation . 

The  improvement  in  diagonal  dominance  of  the  truncated 
matrices  can  be  seen  in  appendix  C.  Not  only  are  the  row  sum  norms 
less  for  Optimal  Weighting  Functions  than  for  unmodified 
Biorthogonal  Weighting  functions,  but  they  are  decreasing  with  the 
row  index,  and  they  are  less  subject  to  the  effects  of  truncation. 

As  an  example  of  the  increased  stability  in  the  early 
coefficients,  the  first  two  coefficients  for  all  the  orders  of 
truncation  shown  in  Appendix  D  for  0__  *  l-2ra,  o__  a  0  are 

as  A  Z 
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N®5 

-0 . 11902E- 1 

0. 11986E-1 

N=iO 

-0. 16985E-1 

0. 14775E-1 

N=  20 

-0 . 16622E-1 

0. 14S78E-1 

N=50 

-0. 16582E-1 

0 . 145S6E- 1 

Ns  100 

-0 . 16566E-1 

0. 14S47E-1 

O . 34817E- 3  -O . 16901E- 3 
0.47215E-4  -0.2524OE-3 
0 . 66895E-4  -0.24615E-3 
0.69037E-4  -0.24546E-3 
0 . 69904E-4  -0.2451BE-3 


for  Unmodified  Bi orthogonal  Weighting  Functions,  and 


C 

1 

C 

a 

N=5 

-0. 16470E-1 

0 . 14509E-1 

0.75161E-3 

-0.2425 IE- 3 

N=  10 

-0. 16588E-1 

0. 14566E-1 

0.6B0S7E-4 

-O . 24554E- 3 

Ns  20 

-0. 16572E-1 

0. 14549E-1 

0 . 69270E-4 

-0.24531E-3 

Ns  SO 

-0. 165SBE-1 

0 . 14543E- 1 

0 . 70286E-4 

-0. 24505E-3 

N=  lOO 

-0. 165S7E-1 

0 . 14S42E-1 

0 . 70409E-4 

-0 . 24S02E—  3 

for  Optimal  Weighting  functions .  The  corresponding  coefficients 


for  the  incompatible  distribution  a 


O,  tr  *  r,  which  presents 


l  more 

severe  test 

of  convergence 

and  stability,  are 

C 

c 

l 

a 

Ns  5 

-0 . 10644E+0 

0.S7713E-1 

— 0 . 6B485E— 2  -0.14244E-2 

N*  lO 

-0. 31706E-1 

0 . 26662E- 1 

-0 . 23732E— 2  -0.20218E-3 

Ns  20 

-0.1142BE-1 

0. 15645E-1 

-0 . 12832E-2  O . 14B38E—  3 

N»50 

0 . 83730E- 2 

0 . 48957E-2 

-0 . 22527E- 3  0.49061E-3 

N-100 

0 . 5486SE-1 

-0 , 20350E- 1 

0 . 22660E-2  0.12874E-2 

for  Unmodified  biorthogonal  weighting  functions  and 


N«5  0.27844E-1  -0.64030E-2 
N* 10  0.29602E-1  -0.72S16E-2 
N«20  0.310B4E-1  -0.79308E-2 
N>50  0 . 32631E-1  -0.B6280E-2 
N«100  0.33316E-1  -0.902S1E-2 


0 . 63268E-3  0.B2867E-3 
0 . 770B1E- 3  0 . 8S079E- 3 
0 . 88376E- 3  0.8749SB-3 
0.9992SE-3  0.902SOE-3 
0 . 1064BE-2  0.9188OE-3 


for  optimal  weighting  function*.  The  increase  in  stability  for  the 
smooth  first  distribution  is  marked,  and  for  the  incompatible  case 


o.w.F.  coefficients  are  still  reasonably  stable,  whereas  the 
u.B.w.F.  coefficients  lose  all  stability. 

The  third  advantage  can  be  seen  in  the  improvement  in 
accuracy  of  the  summed  expansions  tested  against  the  prescribed 
stresses  on  z>0.  Although  the  difference  is  only  slight  for  the 
well-behaved  distribution  l-2r3,  U.B.W.F.  completely  fail  to 
converge  to  the  incompatible  shear  stress,  whereas  the  o.w.F. 
produce  reasonably  good  results  when  the  Cesaro  sums  are 
calculated  rather  than  partial  sums,  as  shown  by  the  graphs  in 
appendix  F . 
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Appendix  A 
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<♦»  ♦«■'>  <"*■» 


«lmk»Vk„>VV  [r-ih-il-<1-r,>x«‘»J.<V'i<v 

1  in  n 

-4k.k»<k.Vk.,Vk»>-k»Vk..,Vk»,l 

-  » 1 kX<  V , <  V J .  < k„> 4k„J .  <  v J ,  <  V  >  1 

(♦n  ♦«  >  <m*n> 

<k.> Ji <  V [(>  »-X  »)ll <1-'>  <ki4ki> J, <k«>J, <k„> 

1  ro  n  ' 

-WV.'VVV^nVVVV’l 

4  7r~a'-T~  .).l2k^n<k.J.<k.)Ji<V-knJi<k.>J.<kr.))l| 

m  n  * 

<*  (ai  .  (4i  v  4  .  % 

♦n  *«  >  «*'“> 

1  n  n 

*  <0^7y7|<k4,',k«k»Ck-',.(k«,J.<kn>-kn'.<k->J.<kn>,)l 


WO 


-  W>«' 


f‘«wj;<v 
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wo 
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lx  »,« 
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WO 
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wo 
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Appendix  B 


This  appendix  lists  explicit  expressions  for  the  right-hand  sides 
dm  corresponding  to  the  six  special  cases  of  section  5,  obtained 
from  optimal  weighting  functions . 
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Appendix  C 


Tl\is  appendix  lists  the  row  sum  norms 

D*m»l  /  l*»J 

n 

for  various  orders  of  truncation  (by  the  order  of  truncation  we 
mean  the  number  of  pairs  of  eigenvalues  used  in  truncating  the 
matrix.  Thus  N= 10  means  that  a  20*20  matrix  has  been  inverted).  It 
should  be  noted  that  although  the  norms  for  the  optimal  weighting 
functions  are  not  less  than  one,  they  decrease  with  n,  (ignoring 
the  effects  of  truncation  for  n  at  either  end  of  the  range)  unlike 
those  for  the  unmodified  biorthogonal  weighting  functions,  and  the 
improvement  this  affords  is  demonstrated  by  the  results  in  the 
next  appendix. 

Biorthogonal  Weighting  Optimal  Weighting 
Functions  Functions 


n 

Ns  5 

N  s  5 

1 

3.60230 

1.71249 

2 

3. 3624B 

2.11606 

3 

3.51319 

1 . 99478 

4 

3 . 295  36 

1 . 83499 

5 

2.16175 

1 . 68600 

z 

II 

H 

o 

O 

H 

II 

Z 

1 

5.61611 

2. 34593 

2 

4.11965 

2.88797 

3 

4 . 22286 

2.77801 

4 

4.46599 

2.61097 

5 

4.68817 

2.44794 

6 

4.86089 

2. 29953 

7 

4.95964 

2.16658 

8 

4.90947 

2.04769 

9 

4.50296 

1 . 94108 

10 

3. 12498 

1.84507 

N  a  20 


K  *  20 


1 

2 

3 

4 

5 

6 
7 
B 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 


9.59345 
5.47667 
4.91307 
4.90790 
5 . 07706 
5 . 30875 
5 .56248 
5.81666 
6 . 06628 
6.30916 
6 .54364 
6.76732 
6 . 97645 
7.16500 
7 . 32200 
7.42656 
7.43424 
7 . 24075 
6 . 57820 
4.77745 


3.00200 
3.69238 
3.60346 
3.44068 
3.27596 
3.12280 
2.98317 
2.85629 
2.74077 
2.63518 
2 . 5  3827 
2.44895 
2. 36631 
2.28958 
2. 21808 
2.15128 
2.08868 
2.02987 
1.97450 
1 . 92225 


N  *  50 


N  *  50 


2 

9.55000 

4 

6.15013 

6 

5.87586 

8 

6.15164 

10 

6.58065 

12 

7.05966 

14 

7.55387 

16 

8.04895 

18 

8.53647 

20 

9.01408 

22 

9.48194 

24 

9.94036 

26 

10.38957 

28 

10 . 82965 

30 

11 . 26041 

32 

11.68129 

34 

12.09109 

36 

12.48757 

38 

12.86643 

40 

13.21917 

42 

13.52738 

44 

13.74547 

46 

13.73585 

48 

12. 95004 

50 

8.77780 

4 . 77807 
4.56876 
4.25654 
3.98821 
3.76178 
3.56802 
3.39950 
3,25085 
3.11817 
2 . 99860 
2.88994 
2.79052 
2.69903 
2.61440 
2.53578 
2.46245 
2.  39383 
2.32943 
2.26881 
2.21161 
2.15753 
2.10628 
2.05763 
2.01136 
1.96730 


30 


N 


100 


N  »  100 


4 

8. 23533 

5.43116 

8 

6.70562 

4.85933 

12 

7.30349 

4.43930 

16 

8.18078 

4.11988 

20 

9.10390 

3.86388 

24 

10.02212 

3.65091 

28 

10.92245 

3.46894 

32 

11.80173 

3.31037 

36 

12.65852 

3.17007 

40 

13.49438 

3.04447 

44 

14. 31146 

2.93093 

48 

15 . 11152 

2.82748 

52 

15 . 89599 

2.73260 

56 

16 . 66597 

2.64509 

60 

17.42227 

2.56399 

64 

18 . 16532 

2.48851 

68 

18.89513 

2.41800 

72 

19.61101 

2. 35193 

76 

20.31100 

2.28983 

80 

20.99060 

2.23130 

84 

21.63895 

2.17603 

88 

22.22654 

2.12370 

92 

22.64989 

2 . 07408 

96 

22. 32804 

2.02692 

lOO 

14.21699 

1 . 98204 

Appendix  D 


Coefficients  and  Summed  Expansions 


convergence  to  om _ 
sz 


i-*ra,  orz  *  o 


Biorthogonal  Weighting 
Functions 


Optimal  Weighting 
Functions 


n 


N  =  5 


N  =  5 


1 

2 

3 

4 

5 


-0 . 11902E- 1 
0. 34817E- 3 
0.88271E-4 
0 . 22726E-4 
-0. 10792E-4 


0. 11986E-1 
0. 16901E-3 
0 . 22602E-4 
0. 23230E-4 
0. 39198E-4 


-0. 16468E-1 
0.75162E-4 
0. 32191E— 4 
0 . 76576E-5 
0.21059E-S 


O. 14509E-1 
-O. 2425 IE- 3 
-0. 10014E-4 
0. 13609E-5 
0. 11227E-5 


r  Normal  stress  shear  stress 


°2Z 

UBWF 

OWF 

arz 

UBWF 

OWF 

0.0 

1.00 

0.8283 

1.0243 

0.00 

0.0000 

0  .  OOOO 

0.1 

0.98 

0.8834 

0 . 9830 

0.00 

-0.1076 

0.0183 

0.2 

0.92 

0.8936 

0 . 9070 

0.00 

-0.0177 

-0.0036 

0.3 

0.82 

0.7422 

0 . 8277 

0.00 

0.0817 

-0.0121 

0.4 

0.68 

0.5895 

0 . 6882 

0.00 

0.0376 

0.0130 

0.5 

0.50 

0.4723 

0.4841 

0.00 

0.0017 

0 . 0087 

0.6 

0.28 

0 . 2426 

0.2781 

0.00 

0.0414 

-0.0180 

0.7 

0.02 

-0.0451 

0 . 0386 

0.00 

0.0621 

0.0007 

0.8 

-0.28 

-0.2866 

-0.2892 

0.00 

0.0896 

0 . 0228 

0.9 

-0.62 

-0.5559 

-0.6262 

0.00 

0.0965 

-0.0051 

1.0 

-1.00 

-0.7252 

-1.0034 

0.00 

0  .  OOOO 

O . OOOO 

NOTE 

ubwf  -  unmodified  Biorthogonal  Weighting  Functions 
OWF  -  optimal  Weighting  Functions 


Bi orthogonal  Weighting 
Functions 


optimal  Weighting 
Functions 


N  =  10 


N  =  10 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


-0. 16985E-1 
0.47215E-4 
0. 27679E-4 
0.6S262E-5 
O. 17761E-5 
0 . 59203E-6 
0. 27527E-6 
0 . 19793E-6 
0. 12656E-6 
-O . 96880E- 7 


0 . 14775E-1 
-0.25240E-3 
-0. 13361E-4 
0 . 67046E-7 
0 . 5275SE-6 
0 . 27933E-6 
0. 13833E-6 
0 . 11598E-6 
0 . 19600E-6 
0. 13266E-7 


-0 . 16588E-1 
0. 6B058E-4 
0. 30973E-4 
0 . 73391E- 5 
0. 20022E-5 
0. 62190E-6 
0 , 21096E-6 
0.74706E-7 
0. 26008E-7 
O. 78339E-8 


0 . 145S6E- 1 
-0.245S4E-3 
-O. 11049E-4 
O. 9S283E-6 
0. 93628E-6 
0 . 49396E- 6 
O . 24802E- 6 
0.12711E-6 
0.67366E-7 
O. 36903E- 7 


r  Normal  Stress 


Shear  stress 


°2Z 

UBWF 

0.0 

1.00 

1.0173 

O  .  1 

0 . 98 

O . 9825 

0.2 

0.92 

0.9298 

0.3 

0.82 

0.8237 

0.4 

0.68 

0.6874 

0.5 

0.50 

0.5047 

0.6 

0.28 

0.2839 

0.7 

0.02 

0.0241 

0.8 

-0.28 

-0.2805 

0.9 

-0.62 

-0.6254 

1.0 

-1.00 

-1.0257 

OWF 

a 

UBWF 

rz 

0. 9794 

0.00 

0 . OOOO 

0 . 9880 

0 . 00 

0.0032 

0. 9135 

0.00 

-0.0032 

0.8264 

0.00 

-0.0013 

0 . 6744 

0.00 

-0.0027 

0.5056 

0.00 

-0.0038 

0 . 2754 

0.00 

-0.0043 

0.0239 

0.00 

-0.0045 

0.2817 

0.00 

-0.0083 

0.6213 

0.00 

-0.0049 

0.9990 

0.00 

0.0000 

OWF 

0.0000 
-0.0038 
0.0012 
0.0002 
-0.0019 
0.0029 
-0.0050 
0.0058 
-0.0081 
0.0072 
O . OOOO 
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B  ior thogonal 

Weighting 

Optimal  weighting 

Functions 

Functions 

N  = 

20 

N  = 

20 

1 

-0.16622E-1 

0. 14578E-1 

-0 . 16572E-1 

0.14549E 

2 

0 . 66895E-4 

-0.24615E-3 

0 . 69270E-4 

-0.24531E 

3 

0. 30920E-4 

-0. 11202E-4 

0. 3125 IE-4 

-0.10945E 

4 

0. 73657E-5 

0 . 91154E- 6 

0 . 74367E- 5 

0. 99956E 

5 

0 . 20277E-5 

0 . 92480E-6 

0 . 20454E- 5 

0 . 96002E 

6 

0 . 64032E-6 

0 . 49145E- 6 

0 . 64404E-6 

0 . 5  0723E- 

7 

0 . 22392E-6 

0. 24839E-6 

0. 22349E-6 

0. 25600E- 

8 

0. 84011E-7 

0. 12839E-6 

0 . B2323E-7 

0. 13220E 

9 

O. 32906E-7 

0. 68881E-7 

0. 30901E-7 

0.70754E- 

10 

0. 13125E-7 

0. 38454E- 7 

0. 11116E-7 

O. 39245E 

11 

0 . 52022E-8 

0. 22364E-7 

O. 3294 9E-8 

0. 22481E 

12 

O . 19950E-8 

0. 13609E-7 

0 . 23007E-9 

0. 13241E 

13 

0.69736E-9 

0 . 87645E- 8 

-0 . 88613E-9 

0.79837E 

14 

0. 12389E-9 

0 . 60958E-8 

-0 . 11986E-8 

0.49067E 

15 

-0. 30167E- 9 

0 . 46801E- 8 

-0. 11890E-8 

O. 306 10E 

16 

-0. 97063E-9 

0. 39408E- 8 

-0 . 10627E-8 

0.19300E 

17 

-0 . 22738E-B 

0. 32481E-8 

-0. 90673E-9 

0.12244E 

18 

-0 . 43301E-8 

0. 12425E-8 

-0.75622E-9 

O . 77744E 

19 

-0.43209E-8 

-0.42482E-8 

-0.62369E-9 

0. 49075E 

20 

0. 19237E-8 

-0.76306E-9 

-0 . 51189E- 9 

O. 30513E 

r 

Normal  Stress 

Shear  stress 

°ZZ 

UBWF 

OWF 

°rz 

UBWF 

OWF 

0.0 

1 . 00 

0.9985 

0. 9937 

0.00 

0.0000 

0.0000 

0.1 

0.98 

0.9803 

0 . 9785 

0.00 

0.0006 

0.0012 

0.2 

O.  92 

O. 9204 

0 . 9188 

0 . 00 

0.0001 

0 . 0008 

O.  3 

0.82 

0.8206 

0.8188 

0.00 

-0.0003 

0.0006 

0.4 

0.68 

0.6807 

0.6788 

0.00 

-0.0005 

0 . 0005 

0.5 

0.50 

0.5008 

0 . 4986 

0 . 00 

-0.0006 

0.0003 

0.6 

0.28 

0.2809 

0.2784 

0.00 

-0.0006 

O .  OOOl 

0.7 

0.02 

0.0208 

0.0182 

0.00 

-0.0006 

-0.0002 

0.8 

-0 . 28 

-0. 2798 

-0 . 2820 

0.00 

-0.0006 

-0.0007 

0.9 

-0.62 

-0 . 6210 

-0.6218 

0 . 07 

-0 . 0010 

-0.0014 

1.0 

-1.00 

-1.0039 

-0.9995 

0.0  0 

0 . OOOO 

O .  OOOO 
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Biorthogonal  Weighting 
Functions 


Optimal  Weighting 
Functions 


N  =  50  N  S  50 


2 

0 . 69037E-4 

-0. 24546E-3 

0 . 70286E-4 

-0.24505E-3 

4 

0. 74558E-5 

0. 99827E-6 

0 . 75074E-5 

O. 10471E-5 

6 

0.65137E-6 

0 . 5 1019E-6 

0 . 65842E-6 

0.S2040E-6 

Q 

0 . 85310E-7 

0. 13409E-6 

0 . 86S78E-7 

0. 13720E-6 

io 

0. 12510E-7 

0.40352E-7 

0 . 12961E- 7 

0.41542E-7 

12 

0. 95515E-9 

0 . 1 3909E-7 

0 . 11042E-8 

O. 14440E- 7 

14 

-0. 78839E-9 

0.53280E-8 

-0 . 73710E-9 

O. 55911E-8 

16 

-0 . 81470E- 9 

0 . 22066E-8 

-0  . 79879E- 9 

0. 23482E-8 

18 

-0 . 59796E- 9 

0. 96547E-9 

-0.5  9576E- 9 

0. 10466E-8 

20 

-0 . 40622E- 9 

0. 43691E-9 

-0. 40949E-9 

0 . 48577E- 9 

22 

-0. 27090E-9 

0. 19987E-9 

-0 . 27630E-9 

0. 23048E-9 

24 

-0. 18075E-9 

0 . 89619E- 10 

-0 . 18689E- 9 

0 . 10939E- 9 

26 

-0 . 12143E-9 

0. 37238E- 10 

-0. 12773E-9 

0 . 50267E- 9 

28 

-0 . 82238E- 10 

0. 12297E-10 

-0 . 88466E-10 

0. 2O936E-10 

30 

-0.5608  3E- 10 

0.74857E-12 

-0 . 62144E-10 

0 . 63731E-11 

32 

-0. 3B427E- 10 

-0 . 41172E- 11 

-0.44268E-10 

-0 . 70643E-12 

34 

-O. 26417E-10 

-0 . 55884E- 11 

-0. 31959E-10 

-0. 39460E-11 

36 

-0. 18298E-10 

-0.5  3050E-11 

-0. 23364E-10 

-0. 52122E-11 

38 

-0. 1 3067E- 10 

-0 . 41206E- 11 

-0. 17283E-10 

-0.54785E-11 

40 

-0. 10321E-10 

-0 . 26 1 05 E- 11 

-0. 12956E-10 

-C , 52578E-11 

42 

-O. 10188E-10 

-O. 16237E-11 

-0. 97646E-11 

-0 . 48209E- 1 1 

44 

-0 . 12881E- 10 

-O. 34865E- 11 

-0 . 74462E- 11 

-0.4  3096E- 1 1 

46 

-O. 14402E-10 

-0. 14170E-10 

-0.57276E-11 

-0. 37959E-11 

48 

0. 13179E-10 

-0. 24198E-10 

-  0 . 4 44 HE- 11 

-0. 33145E- 1 1 

50 

-O . 5B778E- 11 

O . 5  3884E- 11 

-0. 34692E-11 

-0.28799E-11 

r  Normal  stress  Shear  Stress 


UBWF 

OWF 

UBWF 

OWF 

0.0 

1.00 

1 . 0006 

1.0000 

0.00 

0 . 0000 

0.0000 

0.1 

0.98 

0.9803 

0 . 9800 

0.00 

o . oooo 

-0.0001 

0.2 

0.92 

0.9204 

0.9201 

0.00 

-0.0001 

0.0001 

O.  3 

0.82 

0.8204 

0.8200 

0 . 00 

-0.0001 

-0.0002 

0.4 

0.68 

0.6803 

0.6800 

0.00 

-0.0002 

0.0002 

0.5 

0.50 

0.5003 

0.5000 

0.00 

-0.0002 

-0.0003 

0.6 

0.28 

0.2802 

0.2799 

0.00 

-0.0003 

0 . 0004 

0.7 

0.02 

0.0202 

0.0202 

0.00 

-0.0003 

-0.0005 

0.8 

-0.28 

-0.2800 

-0.2803 

0.00 

-0.0004 

0 . 0005 

0.9 

-0.62 

-0.6203 

-0.6195 

0 . 00 

-0.0005 

-0.0006 

1.0 

-1.00 

-1.0015 

-0.9999 

0.00 

0 . 0000 

0.0000 
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Biorthogonal  weighting 
Functions 


N  =  100 


4 

0 . 74939E-5 

0. 10328E-5 

8 

0 . 86652E-7 

0. 1364SE-6 

12 

0. 11147E-8 

0. 14343E-7 

16 

-0.78652E-9 

0. 23318E-8 

20 

- 0  .  5025 OE- 9 

0 . 48385E- 9 

24 

-0 . 18301E- 9 

0. 11034E-9 

28 

-0 . 86235E- 10 

0. 22270E-10 

32 

- 0 . 42922E- 10 

0.45  789E- 12 

36 

-O. 22513E-10 

-0 . 4 300 IE- 11 

40 

-0. 12360E-10 

-0.45703E-11 

44 

-0.70510E-11 

-0. 38008E- 11 

48 

~0 . 41486E- 11 

-0. 29431E-11 

52 

-0. 24975E-11 

-0. 2229SE-11 

56 

-0. 15233E-11 

-0. 16856E- 11 

60 

-0. 92723E-12 

-0. 12842E-11 

64 

-0 . 547S8B— 12 

-0. 99108E-12 

68 

-0. 29289E-12 

-0. 77628E-12 

72 

- 0 . 10869E- 12 

-0.6155  3E- 12 

76 

O. 39932E- 1 3 

-0 . 48769E-12 

80 

O. 17718E-12 

-0. 36936E-12 

84 

0. 31717E-12 

-0 . 22381E-12 

88 

O . 4  3820E- 12 

0 . 22179E- 12 

92 

0. 32371E-12 

0 . 47258E- 12 

96 

-O . 79584E- 12 

0 . 27566E- 12 

lOO 

-0.77300E-13 

0 . 15096E- 12 

Optimal  Weighting 
Functions 


N  =  lOO 

0 . 75161E-5  0 . 10S29E-5 

0 . 87440E- 7  0 . 1 3762E-6 

0.12125E-B  0.14593E-7 

- 0 . 76611E- 9  0 . 24035E-8 

-0 . 39691E- 9  0 . 5 1042E- 9 

- 0 . 18125E- 9  O . 1220 IE- 9 

-0 . 85654E-10  0.280S3E-10 

-0 . 4 275 5 E- 10  0.36001E-11 

-0.22503E-10  - 0 . 24601E- 1 1 
- 0 . 12412E- 10  -0.34212E-11 
- 0 . 71294E- 1 1  - 0 . 30400E- 1 1 
- 0 . 42400E- 1 1  -0.24108E-11 
-0.25980E-11  -0.18365E-11 
— 0. 16333E-11  -O. 13798E-11 
- O . 10496E- 11  -0 . 10343E-11 
-0.68737E-12  -0.77781E-12 
-0.45744E-12  -0.58833E-12 
-0.30857E-12  -0.44819E-12 
- 0 . 21049E- 12  -0.34406E-12 
-0.14488E-12  -0.26617E-12 
-O . 10039E- 12  -0.20750E-12 
-0.69865E-13  -0.16297E-12 
-O . 48708E- 1 3  -0.12890E-12 
-0.33916E-13  -0.10265E-12 
-O . 23500E- 1 3  -0.82266E-13 


r  Normal  stress  shear  Stress 


°zz 

UBWF 

OWF 

axz 

UBWF 

OWF 

0.0 

1.00 

1.0027 

0.9997 

0.00 

0 . 0000 

O . OOOO 

0.1 

0.98 

0 . 9804 

0.9800 

0.00 

-0.0003 

O . OOOO 

0.2 

0 . 92 

0.9204 

0.9200 

0.00 

-0.0002 

0 . OOOO 

0.3 

0.82 

0.8204 

0.8199 

0.00 

-0.0002 

O . OOOO 

0.4 

0.68 

0.6803 

0.6799 

0.00 

-0.0002 

0 . OOOO 

0.5 

0.50 

0.5003 

0.4999 

0.00 

-0.0002 

0 . OOOO 

0.6 

0.28 

0.2803 

0.2799 

0.00 

-0.0002 

0.0001 

0.7 

0.02 

0.0203 

0.0199 

0.00 

-0.0003 

0.0001 

0.8 

-0.28 

-0.2798 

-0.2801 

0.00 

-0.0003 

0.0002 

0.9 

-0.62 

-0.6200 

-0.6201 

0.00 

-0.0002 

0 . 0002 

1.0 

-1.00 

- 1 . 0006 

-1 . 0000 

0.00 

0.0000 

O . OOOO 
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Convergence  -to  a, 


Of  o 


r 


ZZ 


rz 


NOTE 


The  partial  sums  for  this  stress  distribution  are  not 
convergent,  and  the  sums  shown  below  are  Cesaro  sums, 
i.e. 

n 


If  S  is  the  nth  partial  sum,  then  C  m 


-[  s. 
ni=l  1 


Biorthogonal  weighting 
Functions 


Optimal  Weighting 
Functions 


n 


N  =  5 


N  =  5 


1 

2 

3 

4 

5 


-O. 10644E-0 
-0.6B4S5E-2 
-0 . 13217E-2 
-0. 34379E-3 
0. 20885E-3 


0. 677X3E-1 
-0. 14244E- 2 
-0. 7105 BE- 3 
-0 . 44117E- 3 
-0.S1413E-5 


0. 27844E-1 
0.63268E-3 
0.51659E-4 
0. 39435E-5 
-0. 18B70E-5 


-0 . 64030E- 2 
0 . 82B67E- 3 
O. 16962E-3 
0.S2381E-4 
0.20690E-4 


r  Normal  stress  Shear  stress 


°ZZ 

UBWF 

OWF 

°rz 

UBWF 

OWF 

0.0 

0.00 

3.8987 

0.2517 

0.00 

0.0000 

0.0000 

0.1 

0.00 

2.9468 

0.1779 

0.10 

0.6512 

0.0356 

0.2 

0.00 

1.5725 

0.0983 

0.20 

0.2314 

0.1574 

0. 3 

0.00 

1.5576 

0 .1361 

0. 30 

-0.7526 

0.3124 

0.4 

0.00 

2.1319 

0.1394 

0.40 

-1.0633 

0 . 3889 

0.5 

0.00 

1 . 9386 

0.0377 

0.50 

-0.6903 

0.4027 

0.6 

0.00 

1 . 3380 

0.0267 

0.60 

-0.7798 

0.4898 

0.7 

0.00 

1 . 1321 

0 . 0766 

0.70 

-1 . 5204 

0.6656 

0  .  B 

0.00 

0.7844 

-0.1215 

0.80 

-1.5147 

0.7029 

0.9 

0.00 

-1.8012 

-0.2834 

0.90 

-0.2523 

0.4080 

1.0 

0.00 

-8.6863 

0.4987 

1.00 

0.0000 

0.0000 
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Biorthogonal  Weighting 
Functions 


Optimal  weighting 
Functions 


N  *  10 


H  *  10 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


-0. 31706E-1 
-0 . 23732E-2 
-0.40086E-3 
-0. 10276E-3 
-O. 31990E-4 
-0. 99216E-5 
-O. 15152E-5 
O.21502E-5 
0. 27344E-5 
-O.51105E-6 


0. 26662E-1 
-0.20218E-3 
-0. 16414E-3 
-0.68799E-4 
-0. 31147E-4 
-0. 15341E-4 
-0.76746E-5 
-0.28967E-5 
0. 11829E-5 
0. 190S0E-6 


0.29602E-1 
0.770B1E-3 
0.82S72E-4 
0. 14477E-4 
0 . 26161E-5 
0. 25682E-7 
-0 . S2883E-6 
-0. 57641E-6 
-0.50190E-6 
-0.41006E-6 


-0.72516E-2 
0 . 85079E- 3 
0 . 18016E- 3 
0.57098E-4 
0. 23042E-4 
O. 10861E-4 
0.57046E-5 
0. 32458E-S 
0.19643E-5 
O. 12487E-5 


Normal  Stress 


Shear  stress 


a 

UBWF 

OWF 

Z.Z 

» 

0.0 

0.00 

1.1891 

0 . 0428 

0.1 

0.00 

1.0300 

0.1108 

0.2 

0.00 

1.0541 

0.0698 

0.3 

0.00 

O  .  9870 

0.0946 

0.4 

0.00 

0.9331- 

0.0539 

0.5 

0.00 

0.8375 

0.0659 

0.6 

0.00 

0.6878 

0 . 0282 

0.7 

0.00 

0.4734 

0 . 0087 

0.8 

0.00 

0.1128 

0.0003 

0.9 

0.00 

-0.7781 

-0.1916 

1.0 

0.00 

-3.8650 

0.5033 

a 

UBWF 

OWF 

rz 

0.00 

O .  OOOO 

O . OOOO 

0.10 

0.0311 

0.1071 

0.20 

-0.0559 

0.1745 

0.30 

-0.0084 

0 . 2875 

0.40 

-0.0862 

0.3565 

0.50 

-0.0698 

0.4673 

0.60 

-0.1509 

0.5248 

0.70 

-0.1756 

0.6414 

0.80 

-0.2722 

0.6588 

0.90 

-0.0649 

0.7860 

1.00 

O . OOOO 

0 . OOOO 
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1 


Biorthogonal  weighting 
Functions 


optimal  weighting 
Functions 


N  *  20 


N  =  20 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 


-0. 11428E-1 
-0. 12832E-2 
-O. 22399E-3 
-0. 57686E-4 
-0. 18352E-4 
-0.64103E-5 
-0.21509E-5 
-0 . 46193E-6 
0,25 34 IE- 6 
0.S6S04E-6 
0.69655E-6 
0. 74017E-6 
O, 73143E-6 
0, 67518E-6 
0 . S5031E-6 
0. 30160E-6 
-0. 16928E-6 
-0 . 90938E-6 
-0.11247E-5 
0 . 48445E- 6 


O. 15645E-1 
O. 14838E-3 
-0 . 44775E-4 
-0.23268E-4 
-O. 10654E-4 
-O. 50486E-5 
-0.24762E-5 
-O. 12170E-5 
-O . 5555 9E-6 
-O, 17913E-6 
0.58S45E-6 
0.23229E-6 
O. 38204E-6 
O, 53240E-6 
0.69484E-6 
O . 85477E-6 
O, 91311E-6 
O  .  5 1023E-6 
-O, 10165E-S 
-O. 18472E-6 


Normal  Stress 


O. 31084E- 1 
0 . 88376E- 3 
O, 10765E-3 
0 . 23082E— 4 
0, 63527E-5 
0 . 19086E-5 
O . 52231E-6 
0 . 5594 9E-7 
-O. 99179E-7 
-O . 14183E-6 
-0 . 14318E-6 
-0. 13029E-6 
- O . 11 368E-6 
-O. 973B9E-7 
-0.82821E-7 
-0 . 70304E-7 
-0.59745E-7 
-0.50911E-7 
-0.43540E-7 
-O. 373B8E-7 


-O. 79308E-2 
0 , 87495E- 3 
O. 19064E-3 
0.61745E-4 
0.25 371E-4 
O. 12148E-4 
0.64715E-5 
O. 37300E-5 
O . 22845E- 5 
O. 14685E-5 
O. 98186E-6 
0.67836E-6 
0 . 48185E-6 
O. 35048E-6 
O. 26023E-6 
0.1967JE-6 
0.15111E-6 
0.11772E-6 
0.92877E-7 
0.74119E-7 


Shear  stress 


a  UBWF  OWF 

zz 


rz 


UBWF 


OWF 


0.0 

0.00 

-0.0399 

0.0326 

0.1 

0.00 

0.6073 

0.0500 

0.2 

0.00 

0.6601 

0.0487 

0. 3 

0.00 

0.6704 

0.045  3 

0.4 

0.00 

0.6528 

0.0403 

0.5 

0.00 

0.6036 

0.0340 

0.6 

0.00 

0.5086 

0.0262 

0.7 

0.00 

0.3382 

0 . 0170 

0.8 

0.00 

0.0342 

0.0053 

0.9 

0.00 

-0.5280 

-0.0089 

1.0 

0.00 

-2.4519 

0.5068 

0.00 

0 . OOOO 

0.0000 

0.10 

0 . 2237 

0.0883 

0.20 

0 . 2075 

0.1830 

0.  30 

0.2054 

0.2760 

0.40 

0 . 2007 

0 . 3677 

0.50 

0. 1868 

0.4577 

0.60 

0.1604 

0.5452 

0.70 

0.1226 

0.6291 

0.80 

0.0893 

0.7049 

0.90 

0.1204 

0.7490 

1.00 

0 . 0000 

O . 0000 

s. 


-  39  - 


I 


i 

+ 


Bi orthogonal  Weighting 
Function* 


optimal  weighting 
Function* 


N  =  50  N  *  50 


2 

-0.22527E-3 

0.49061E-3 

0. 99925E-3 

0. 902SOE-3 

4 

-0.13663E-4 

0 . 19719E-4 

0 . 31656E-4 

0 . 66934E-4 

6 

-0. 12284E-5 

0 . 42449E-5 

0. 37754E-5 

0.13609E-4 

8 

0. 16323E-7 

0. 16176E-5 

0. 68948E-6 

0 . 42925E-5 

10 

0 . 15705E-6 

0 . 78827E-6 

0. 13239E-6 

0. 17304E-S 

12 

0. 14700E-6 

0.44161E-6 

0. 80056E-8 

0 . 81676E— 6 

14 

0. 11961E-6 

0. 27186E-6 

-0. 19961E-7 

0 . 43045E-6 

16 

0. 9S898E-7 

0. 17971E-6 

-0. 23538E-7 

0.24612E-6 

18 

0.77936E-7 

0. 1259SE-6 

-0 . 2099SE-7 

0. 14984E-6 

20 

0.64636E-7 

0. 92907E-7 

-0. 1736BE-7 

0. 95888E-7 

22 

0.54708E-7 

0.71878E-7 

-0. 14029E-7 

0 . 63900E-7 

24 

0.47131E-7 

0 . 5B240E-7 

-0. 11274E-7 

0.44040E-7 

26 

0.41151E-7 

0 . 4942SE-7 

-0. 90829E-8 

0. 31226E-7 

28 

0. 36190E-7 

0.43949E-7 

-0.73614E-8 

0. 22684E-7 

30 

0. 31753E-7 

0 . 40948E-7 

-0.60103E-8 

O. 16828E-7 

32 

0. 27338E-7 

0. 39922E-7 

-0.49457E-8 

O . 12715E-7 

34 

0.22319E-7 

0.40575E-7 

-0 . 41015E-8 

0. 97631E-8 

36 

0.1S778E-7 

0 . 42660E-7 

-0. 34273E-8 

0 . 76050E-8 

38 

0.62169E-8 

0 . 4S703E-7 

-0. 28846E-8 

0.60004E-8 

40 

-0.88936E-8 

0.48255E-7 

-0 . 24444E-8 

0 . 47893E- 8 

42 

-0. 33468E-7 

0 . 4S68SE-7 

-0.20846E-8 

0. 38628E-8 

44 

-0.70257E-7 

0. 23722E-7 

-0. 178BSE-8 

0. 31452E-8 

46 

-O. 95703E-7 

-O.S2967E-7 

-0. 15431E-8 

0.25B32E-8 

48 

0.5051BE-7 

-0 . 13265E-6 

-0. 13384E-8 

0.21385E-8 

50 

-0.2726SE-7 

0. 25884E-7 

-0. 11665E-8 

0 . 17834E-B 

r  Normal  stress  Shear  stress 


°*Z 

UBVF 

OWF 

°rz 

UBWF 

OWF 

0.0 

0.00 

2.0988 

0 . 0225 

0.00 

0 . 0000 

0.0000 

0.1 

0.00 

0.2693 

0.0344 

0.10 

0 . 3377 

0.0965 

0.2 

0.00 

0.4966 

0.0300 

0.20 

-0.1055 

0.1888 

0.3 

0.00 

0.4078 

0.0321 

0.30 

0 . 3475 

0.2857 

0.4 

0.00 

0. 3462 

0.0267 

0.40 

0 . 0844 

0 . 3774 

0.5 

0.00 

0.4080 

0 . 0274 

0.50 

0. 3589 

0.4738 

0.6 

0.00 

0.2145 

0.0208 

0.60 

0.2812 

0.5625 

0.7 

0.00 

0.2317 

0.0165 

0.70 

0 . 3425 

0.6595 

0.8 

0.00 

0.0154 

0.0089 

0.80 

0.4142 

0.7390 

0.9 

0.00 

-0.2509 

-0.0397 

0.90 

0 . 4407 

0.8420 

1.0 

0.00 

-0.8829 

0.S091 

1.00 

0 . 0000 

0 . OOOO 
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<D  * 


Biorthogonal  weighting 
Functions 


optimal  weighting 
Functions 


N  »  100 


N  *  100 


12 

16 

20 

24 

28 

32 

36 

40 

44 

48 

52 

56 

60 

64 

68 

72 

76 

80 

84 

88 

92 

96 

lOO 


0. 95239E-4 
0. 37014E-5 
0.52896E-6 
0.13384E-6 
0. 47464E-7 
0 . 21179E-7 
O. 11166E-7 
O . 66574E-8 
0 .43386E-8 
O. 30025E-8 
O. 21476E-8 
O . 15405E-8 
0.105878-8 
0.62831E-9 
O. 19478E-9 
-0.29348E-9 
-O . 89764E-9 
-0.17021E-8 
-0.28297E-8 
-0.44480E-8 
-0.66801E-8 
-0.891 39E-8 
-0 , 56161E-8 
0.2345 IE- 7 
O. 24496E-8 


O . 11908E- 3 
0.83793E-5 
O. 16646E-5 
O. 52082E-6 
0.21090E-6 
0, 10102E-6 
0.54545E-7 
O. 32276E-7 
O. 20572E-7 
O. 13973E-7 
O. 10047E-7 
0.76165E-8 
0. 60709E-8 
0.50753E-8 
0.44357E-8 
0.40310E-B 
O. 37738E-8 
O. 35772E-B 
O. 33077E-8 
O . 26844E-8 
0.10234E-8 
-O. 33647E-8 
-O. 13141E-7 
-0.53547E-8 
-0.46558E-8 


O. 36441E-4 
0. 10350E-5 
O . 83162E-7 
0.21734E-8 
-0.61255E-8 
-0 . 55460E-8 
-0.41231E-8 
-0.29718E-8 
-0 . 215  38E-8 
-O. 158S4E-8 
-0 . 11880E-8 
-O . 90606E-9 
-0.70246E-9 
-0.S5285E-9 
-0 .44107E-9 
-O. 35626E-9 
-0.29097E-9 
- 0 . 24007E- 9 
-0.19989E-9 
-0.16784E-9 
-O . 1420 IE- 9 
-0.12100E-9 
-0 . 10377E-9 
-0.89S27E-10 
-0.77664E-10 


0.69952E-4 

0.46246E-5 

O. 9004 3E- 6 

0.27684E-6 

O . 10983E-6 

0.51 30OE- 7 

0.26B39E-7 

0.15265E-7 

O . 92562E-8 

0.59O49B-8 

O.  3925 3E- 8 

0.26999E-8 

O . 19111E-8 

O. 13863E-8 

O.10272E-8 

0.77525E-9 

0.S9472E-9 

0.462878-9 

O. 36493E-9 

0.29107E-9 

0.23462E-9 

0.19092E-9 

0.15672E-9 

0.12968E-9 

O. 10810E-9 


Normal  stress 


shear  stress 


0 

UBWF 

OWF 

ZZ 

0.0 

0.00 

-3.2079 

0.0273 

0.1 

0.00 

-0.4719 

0 . 0226 

0.2 

0.00 

-0.3235 

0.0222 

O.  3 

0.00 

-0.2290 

0.0216 

0.4 

0.00 

-0.1588 

0.0208 

0.5 

0,00 

-0.1109 

0.0198 

0,6 

0.00 

-0.0881 

0.0184 

0.7 

0.00 

-0.0808 

0.0160 

0.8 

0.00 

-0.0272 

0.0113 

0.9 

0.00 

0.2336 

-0.0020 

1.0 

0.00 

2.7601 

0.5116 

0 

UBWF 

OWF 

r* 

0.00 

0 . 0000 

0.0000 

0.10 

0.4175 

0.0967 

0.20 

0.4858 

0.1926 

0.  30 

0 . 5616 

0.2884 

0.40 

0.6277 

0 . 3839 

0.50 

0.6841 

0 . 4790 

0.60 

0.7432 

0.5734 

0.70 

0.8352 

0.6667 

0.80 

1.0034 

0.7579 

0.90 

1.1914 

0.8412 

1.00 

0.0000 

0 . 0000 
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APPENDIX  E 


Graphical  Results 

The  following  are  a  selection  of  results  obtained  for  the  six 
special  cases  discussed  in  section  S,  all  obtained  by  using 
optimal  Weighting  Functions  and  truncating  the  infinite  matrix 
using  100  eigenvalues. 

Perhaps  the  most  striking  feature  is  the  improvement  in 
convergence  when  Cesaro  sums  are  used  rather  than  partial  sums  in 
those  cases  where  the  coefficients  decay  too  slowly  for  the 
expansions  to  convergence  normally.  As  with  ordinary  Fourier 
series  these  discontinuous  cases  show  Gibbs  phenomena  in  the 
neighbourhood  of  the  jumps,  as  pointed  out  by  Joseph  and 
sturges  (1978)  for  the  semi- infinite  strip. 
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NON-CANONICAL  PROBLEfl 


Convergence  to  the  data  For  z=0 


NON-CANONICAL  PROBLEfl 
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NON-CANONICAL  PRQBLEH 


Convergence  to  the  data  For  z=0 
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Convergence  to  the  data  For  z=0 


NON-CANONICAL  PROBLEM 


Partial  Suns*  100  Terras 


Convergence  to  the  data  For  z=0 


Discontinuous  Shear  Stress  Distribution 


NON-CANONICAL  PROBLEM 


APPENDIX  F 


Note*  on  the  Computations 

All  the  computation®  discussed  in  thia  r®port  ware  carried 
out  on  th®  Oxford  univaraity  Engineering  science  department • a  VAX 
11/780  machin®. 

The  programmes  to  calculat®  th®  ®ig®nvalu®a  and  all  the 

Beaa®l  functiona  required  for  aumming  axpanaiona  were  calculated 

uaing  a  alightly  modified  version  of  the  BRL  Beaael  function 

aubroutin®.  The  programn®  waa  rewritten  in  DOUBLE  complex 

(COMPLEX* 16 )  arithmetic,  and  waa  simplified  alightly  ao  that  only 

the  Baaael  functiona  J  (a)  and  J  (a)  would  b®  calculated  for  each 

•  l 

call  to  th®  aubroutin®.  Th®  ®ig®nvalu®a  could  b®  calculated  for 
any  value  of  Poiaaon'a  ratio  v,  but  all  th®  results  in  thia  report 
hav®  ua®d  th®  value  v  «  0.3.  Th®  eigenvalues,  their  Beaael 
functiona  ( X^)  and  (VJ ,  and  B®aa®l  functiona  of  the  form 
<?a(Xmr)  for  various  intermediate  values  of  rc[0,l],  were 
calculated  in  advance  and  stored  on  disk. 

The  cylinder  eigenvalues  were  calculated  using  a  simple 
Newton  iteration  technique  which  was  found  to  produce  satisfactory 
convergence  to  values  which  agreed  to  virtually  full  double 
precision  with  those  calculated  at  BEL.  The  programme  for 
calculating  the  coefficients  for  the  non-canonieal  stress  problem 
was  built  around  the  NAS  library  routine  F04ADF  which  solves 
complex  systems  of  linear  equations  (with  multiple  right-hand 
sides  if  required)  using  the  Crout  factorisation  method.  The 
subroutines  to  set  up  the  infinite  matrix  and  the  right-hand  sides 


-  93 


were  coded  to  test  both  unmodified  biorthogonal  weighting 
functions  and  optimal  weighting  functions.  The  matrix  was  checked 
for  diagonal  dominance  and  the  equations  were  inverted.  Having 
obtained  the  coefficients,  the  run  could  be  terminated  if  desired, 
otherwise  the  eigenfunction  expansions  could  be  summed  either  for 
a  few  points  in  the  range  (0,1)  to  test  the  convergence  to  the 
prescribed  data  or  over  a  large  number  of  points  for  various 
numbers  of  terms  and  for  increasing  values  of  a  for  use  in  a 
graphics  program. 
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